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The amplitude and phase of the groundwave are calculated for oblique propagation 
across a flat lying coastline. The land and sea are assumed to be smooth and homogeneous. 
Attention is focused on the effects which take place near the coastline. It is shown that the 
reflected wave depends critically on the angle of incidence, O , while the transmitted wave has 
only a weak dependence on 0„. 



1. Introduction 

In certain radio navigational systems it is im- 
portant to estimate the influence of inhomogeneities 
of the earth's surface. A striking example corre- 
sponds to the situation when the transmission is 
over mixed land/sea paths. Here, not only the 
conductivity contrast, but also the change of ele- 
vation of the various portions of the path may 
influence the amplitude and the phase of the received 
signal. 

The problem of calculating fields over mixed 
paths has been investigated by Feinberg [1946], 
Mnlington [1949], Clemmow [1953], Bremmer [1954], 
Furutsu [1956], Wait [1956], Wait and Householder 
[1957], Godzinski [1958], Kalinin and Feinberg 
[1958], and others. Surveys of this and related 
work can be found in a recent review article [Wait, 
1963]. 

It is the purpose of the present paper to consider 
propagation across a straight coastline for oblique 
incidence. Special attention is given to the behavior 
of the fields near the coastline itself since most pre- 
vious theories are not valid there. Also, the reflec- 
tion from the coastline and the change of wavefront 
direction are evaluated. In part II, the influence of 
a gradual elevation change between land and sea is 
treated using an approximate condition due orig- 
inally to Feinberg. Finally, in part III, the effect 
of allowing for the nonabruptness of the effective 
conductivity change between land and sea is 
considered. 

2. Formulation 

To simplify the problem at the outset, a flat coast- 
line is considered and the influence of earth curvature 
is neglected. The situation is illustrated in figure 1. 
With respect to a Cartesian coordinate system, the 
xy plane is taken to be the plane surface of the earth 
and the coastline is the y axis. For purposes of dis- 
cussion, the medium to the left (i.e., x<0) is de- 
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scribed as land with surface impedance Z while the 
medium to the right (i.e., z>0) is described as the 
sen with surface impedance 7J '. The transmitter at 
A with coordinates (—x Q , y ) is regarded as a vertical 
electric dipole, of effective height h a , on the surface 
of the land. The receiving antenna, of effective 
height h b , is located at 13 with coordinates (d u 0) 
where di may be positive or negative. For conven- 
ience, the receiving antenna at B is also located on 
the earth's surface and it is assumed that it responds 
only to the vertical electric field component at 2=0. 
In formulating mixed path problems of this kind, 
it is usually assumed that the surface impedance Z 
for the land, and Z' for the sea, hold right up to the 
boundary. Actually, even if the boundary were ab- 
rupt, there would be a violation of the basic premises 
in the use of the surface impedance concept. The 
conditions for the applicability of the impedance 
relations 

E x =-ZH y and E V =ZH X 

are that Z should change slowly in a distance equal 
to the effective wavelength \ e in the lower medium. 
Because \ e is much less than X , the free space wave- 
length, it means that the surface impedance relations 
may be valid at distances from the boundary small 
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Figure 1. Plane view of the mixed path showing the location of 
the dipoles A and B. 
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compared with X . To extend the usefulness of the 
surface impedance concept, it is desirable to insert 
a transition zone such that Z changes gradually to 
Z'. The width of the transition zone is equal to 2e 
as indicated in figure 1 . It may be defined such that 
Z and Z' are essentially constant when x<^ — e and 
a?> + e, respectively. Apart from mathematical con- 
venience, the existence of such a transition zone has 
a practical significance in studying the influence of 
a variable water depth in the vicinity of the shore- 
line. 

If the problem involved only a homogeneous land 
of surface impedance Z, the mutual impedance z m 
between dipoles A and B would be readily calculated. 
Thus, it is convenient to regard the medium to the 
right of the coastline as a perturbation and the 
resulting change of the mutual impedance is denoted 
Az m . From previous work [Wait, 1956], it is known 
that Az m can be expressed in terms of a surface 
integral over the plane z = 0, of the form 



As. 



-Ml> 



-Z)H a t ' H' bt dS, 



(i) 



where H at is the tangential magnetic field of dipole 

A over the reference earth of constant surface 

-» 
impedance Z and where H' bt is the tangential mag- 
netic field of dipole B over the actual inhomogeneous 
earth. In formulating this integral, the current at 
the terminals of dipoles A and B is /. 

Following earlier work [Wait, 1956], it is conven- 
ient to express the tangential field vectors at the 
variable point P(x, y) m the following forms 



H a 



ikh a I 
~~ 2irr 



( 1+ i) 



e-» r F<r,Z)(i r X%), (2) 



and 



^'=Sf ( 1+ Wi) '—FHB,Z,Z>)&xb, (3) 
where 

r=[(x+x o y+(y-yo) 2 ]i, 

and F and F f are slowly varying quantities. i r , i B , 

-> 
and i z are unit vectors in the directions of increasing 
r, R, and z, respectively. 

The attenuation function F is a known quantity 
from the theory of groundwave propagation over a 
flat homogeneous earth of surface impedance Z. The 
explicit form [Wait, 1963] is given by 



F(r, Z)^l-i(jpfe~ v erfc (ip*), 



(4) 



where p 



__ikr(Z\ 2 



When \p\«l,F(r,Z)o± + l 
while, if \p\^ >l,F(r, Z)o±-l/(2p). The function 



F'(R, Z, Z') is not known except to say that it is a 
function of R, Z, and Z' and that it approaches 
unity if R is sufficiently small. 

-> -> 

Using the above representations for H at and H bt} 
(1) for Az m may be written 



Az m =- 



k 2 h a hi 



4tt 2 



J J 

J X=— € J y= - 



rB 



F(r,Z)F'(R,Z,Z') 



X(Z'- 



(Z '- Z K 1+ ^)( 1+ ^) C ° S ^' (5) 

where 8 is the angle subtended by the units vectors 
-» -> 

i f and i R . Here the element of area dS has been re- 
placed by dxdy. Also, it should be noted that the 
surface S is the region x> — e since Z r — Z^O for 
x< — e. 

3. Approximate Solution 

A number of approximations are now made in 
order to facilitate the integrations. First, it is 
assumed that dipole A is in the far field such that the 
incident wave fronts are approximately plane in the 
neighborhood of the boundary under consideration. 
Then the slowly varying attenuation function 
F(r, Z) may be replaced by F(r , Z) which is the 
value appropriate to the origin. Also, it is assumed 
that the attenuation function F'(R, Z, Z') may be 
replaced by unity if quantities of first order only 
are to be retained. 

Within the approximations stated above it is 
evident that 

cos 5= cos (5 -f-5i)=cos § cos hi — sin 8 sin 8 h 

where 

sin 8 1 =S 1 =y /r , 

cos d 1 = C 1 = (l-S 2 ^=x /r , 

x—di 



sin 8 =y/R and cos 5 =- 



Thus 



R 



cos 5-^ Ci-! S» 



In the phase factor exp(—ik(r-\-R)) it is desirable 
to simplify r by retaining only first-order variations 
in x and y. For example, 



and, therefore, 



r 2 =(x+x ) 2 +(y-y ) 2 
^x 2 o +y 2 o +2xx —2yy 
r^r +xCi—ySi. 



Furthermore, in the integrand of (5), 1/r may be 
replaced by l/r , and the factor 1 + 1 /ikr by unity. 
Using the simplifications indicated, it readily fol- 
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lows thai 



Az», 



4tt 2 



n 



F(r 0i Z) 



X 



/j z '-^-^£"t( i+ «k) 

x [(^ 1 ) cri ~l^ i ] e ^ i% * 



k (6) 



By making use of the identity 
d e~ ikR 



bB E 
it readily follows that 
—ikh a h b e 



•ik 



( 1+ w, 



V 



kRj R 



Az™=- 



4tt 2 






-ikr Q f*ca 

F(r , Z) (7J -Z) exp[-ikC l x] 

r J — • 

e ikS i y dy 



<KL 



n 






dy /* 



The infinite integrals with respect to y may be ex- 
pressed in terms of the Hankel function of order zero 
[Campbell and Foster, 1949]. Explicitly 



I 



e ikS lV fiy- 



j - ikR 



-mn^[kc x \x-<h\i (8) 



which is valid because k can be regarded to have a 
vanishingly small but finite negative imaginary part. 
To simplify the final expression and to cast it into 
a dimensionless form, the results are normalized by 
dividing by the mutual impedance z m , defined by 



ikr\Ji a h b 

ZtqT 



1 



ikr o k 



*_) e -^F(h, Z), (9) 

To/ 



where r =r +d 1 C 1 . This expression corresponds to 
the mutual impedance between dipoles A and B, at 
separation distance r , if they were located on a flat 
homogeneous earth of surface impedance Z. Within 
the limitations already imposed, 



2wr 



e- ilcr °e- iW i d F(r Q , Z). 



(10) 



Making use of (7) and (8), it readily follows that 
Z Z _«./».» 



9 ihC 1 d l /» co 



1o 



X 



/-;*<?, 



0./ 



ff^IJK?,!*-*! 



+TkS 2 1 H ( l ) [kC\ 



wx-d^X 



dx. (11) 



Following an integration by parts, 



2 m 4 



/; [«^-*w] 



X«-«i*ff ( y[Jfca|s-*|]<fe, (12) 

where /(a;) = (Z' — Z)/rj . Introducing the dimension- 
less quantities 

a=kCiX, <xi=kC\di ,and a — fcCie, 

it is evident that 



Az, 
z m 



l -=~e { "i (° H^[\ a - ai \]e- ia 

A J -ao 



(13) 



where f(a) = (Z' — Z)/t) is regarded as a function of a. 

When the vertical stratification in the earth may be 

neglected, the surface impedance function can be 

expressed in terms of the conductivities <r g and a' 

and dielectric constants e* and e '. Thus 



j(*)t 



[{Mh l(<r'+U'c*)l (cr f +i€, W )*J 



where 



e =8.85X10- 12 . 

4. Special Case of Abrupt Boundary 

To effect the integration with respect to a requires 
that the function j\a) is specified. The simplest 
case, of course, is the coastline which can be regarded 
as abrupt. For example, one may consider 



where 



/(«)= -Aoe ir/ *u(a) 

u(a) = l for a>0, 
= fora<0. 



Here, A is to be regarded as a constant and defined 
such that 

Vo 

If both media (i.e., land and sea) are sufficiently 
well conducting that displacement currents can be 
neglected, A is a real quantity. In general, however, 
A is complex and is a measure of the contrast in 
surface impedances between the two homogeneous 
regions. 

For the assumed sharp boundary, it. should be 
noted that 



do 



-A e iir/A d(a), 



where 5(a) is the Dirac impulse function at a—0. 
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Therefore 



Az, 

Zm 



-^-J o °°^ 2 >(|a- ai |) e -^ a ]- 



(14) 



It is immediately evident that this expression for 
the impedance or the field possesses a singularity as 
»! tends to zero. For example, when |«iK<0, 



Hp{\ ai \)i 



2 

a—i- (log \ai\ 

7T 



-log 2+0.5773), (15) 



which would indicate that the field has a logarithmic 
singularity at the boundary. As will be indicated 
below, the singularity is spurious since for any finite 
transition between Z and Z r the field no longer 
behaves in this manner very near the coastline. 
However, provided |«i| is not too small compared 
with unity, the assumption of an abrupt or sharp 
boundary is satisfactory. 

Using the identity 

£ [a«^(ff„«(a)Tiff? ) (a))] = e ±i -ff? ) (a), (16) 
Oct 

it is a straightforward task to show that 



AZrr 
Zm 



=^j/ i3T/i l(C 2 -i<xi)HP(-«i) 



+«,£??> (-a,)]e'«i, (17) 



for ofj-^0, while 



for <*!>(). 

For purposes of computation, it is convenient to 
write 



Az, 

Zm 



l =C 1 A g(a 1 )=A C 1 [flfi(«i)+^ flfctei)]' (19) 



where g x and g 2 do not depend on C x or A . Noting 
that H™(x)=Jm(x)—iYmix) . where J m and Y m are 
Bessel functions of order m, it is a simple matter to 
obtain the following explicit formulas in terms of 
real quantities. 

For x=«i>0, 
0i(x)=2 I cos yf+x) e/oW+sin y-f+^j Yo(x)\ 



+ ^[ sin VT +x ) Jo(x) " 

and 



'3x . 

I ~J-+ X 



) F (x)]> 



(20) 



fc(x)=| [cos (|+x) (J Q -Y 1 )+sm (j+x) (Ji+r )l 

+i | [sin (|+x) (Jo-^O-cos (j+x) («7i+F )]> 

(21) 

where the arguments of the Bessel functions are 
x(=«i) throughout. 
For x=ai<0, 

0i(x)=2 [ cos (t +x ) Jo(— x) + sin (t +x ) F ° (_x) J 

+1 [sin (-f-+ x ) t/"o(— X)— cob (x+ x ) F o(-x)l' 

(22) 
and 

<fe(x)=| [cos (l+x) (Jo+rO-ain (f+x) (Jt-Fo)] 

+* | [sin (|+x) (J +Fi)+cos (|+x) (^-Fo)]' 

(23) 

where the arguments of the Bessel functions are 
— x(— ~ «i) throughout. 

Employing the preceding formulas, the real and 
imaginary parts of g{a x ) are calculated for a range 
of a x from —20 to +20. These results are shown 
in figures 2 to 6 for selected values of O - Additional 
curves of this type may be readily calculated from 
the numerical values of g x and g 2 given in the appendix. 

When the dipole at A is transmitting, the fractional 
change of the field at B, resulting from the presence 
of the medium of surface impedance Z', is (Az m fz m ) 
or CiA g(ai). Thus, the curves in figures 2 and 3 
show the nature of the phenomenon near the coast- 
line. As expected, the field exhibits a singular 
behavior near the boundary. For positive values 
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Figure 4. Real part of the function g(ai) for various oblique 
angles of incidence. 

of a (i.e., B located over the sea), the field is smoothly 
varying and for a x greater than 2 or 3, the imaginary 
part of g(ai) is predominant and it increases to 
relatively large values as a x is increased further. 
For negative values of a h (i.e., B located over the 
land) there is strong evidence that reflection takes 
place at the boundary. Presumably, the standing 
wave pattern results from the interference between 
the incident and the reflected waves. It is interest- 
ing to note that at an angle of 0=45° the reflection 
from the coastline appears to be very weak. This 
fact is confirmed by noting in the asymptotic 
approximation for (— c*i)>>1, that 



; gMi 



V2ir(-ai) 



*[( 1+ iF=5 



32(-«,) ■ ' ' •)} < 2 "> 



It is evident that, if C y ? = cos G d =-^ the leading 

terms in the expansion cancel. In this case, g{a x ) 
is varying as (— ai)~ 3/2 and thus the reflected wave 
would be attenuated quite rapidly. This same 
asymptotic expansion shows that the reflection is 
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Figure 6. Imaginary part of g(ai) for various values of O 
when plotted as a function of (a{)%lcos 2 O . 

strong at highly oblique incidence (where d is 
small). 

The asymptotic behavior of the field for large 
positive values of a t is also very interesting. For 
example, if a x >> 1, 

gM-T^Uci-iaj J^~ e^(i+f+ . . .) 

^i L \ iron \ bai / 



-cti *\i — e 



z3tt/4 I 



('-£+•■■)] *> 

On retaining terms of first order only, 

and if r x denotes the distance, measured in the direc- 
tion of propagation, from the coastline to the ob- 
server, geometry indicates that 

Therefore, within a first order, 

^AoCi^M-iAo^Y, (26) 

I which is independent of the angle Ci. To demon- 
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strate how closely this law is followed in the general 
case, Im g{a x ) for positive a x is plotted as a function 
of (a^/Cl using the complete expressions derived 
abov^e. It is quite evident that for larger values of 
a h the curves for all angles up to 70° approach a 
straight line with a slope of (2/ir)K 

5. Coastal Refraction 

An extremely interesting aspect of oblique propa- 
gation across a coastline is the resulting change of 
wave front direction. To simplify the discussion, it 
is assumed that the field, incident on the boundary 
at x=0, has the form 



E = e~ ikc i x g ikS iv 



(27) 



which is appropriate if the transmitter at A is suffi- 
ciently removed to the left in figure 1. Then the 
resultant field has the form 



-B=£ , (l + A ^(7i)^£ , oexp (A ^Ci), 



(28) 



where A g, defined by (19) et seq., is a function of 
x only. 

By definition the phase velocities in the x and y 
directions are given by 



and 



Vx /c= - (l/k)(d/bx) phase E, (29) 

v v /c=(l/k)(d/i>y) phase E, (30) 



where c is the velocity of light. Then it readily 
follows that 



and 



v x /c^C 1 -(l/k)(b/i>x) Im (GAotf), (31) 

vJc&Sl (32) 



The effective direction d eff of phase propagation is 
then obtained from 



tan 6 ett ^(Vy/v x ), 



or 



tan e ff=tan 6 { 



[ 



. 1 d , 



(Aoi7)]- 



The refraction error 86 is given by 86=6 < . f! — 6 n . 
to within a first order of smallness in 86, 



tan (8 +86) = t&n d - 



1 



cos' 



(33) 

(34) 

Tims, 

(35) 



which leads to the simple but important result that 

(36) 



„„ sin 6 d T . , ~, x 



In terms of the variable a l = kC i d u and for an ob- 
server at #=</,, the refraction error may be expressed 
in the form 



86i 



^C^flmAo^)]. 



(37) 



Using the general expression for g{a x ) given by (19), 
it is found that, for «i>0, 



?=§ Im-fv 



<^), 






(38) 



If A is regarded as real 
5(0) 



^ A -f sin (« x +£) [Jo-CKJo+Yt 



)} 



— COS ( (*! + 



t)W. 



+C(J,-r,)] 



}■ 



(39) 



where, as usual, the arguments of the Bessel func- 
tions are a x . 

At reasonably large distances from the boundary, 
where ai>>l, (38) simplifies to 



iS 



1 \2ttcJ 



Re An 



^sin d^ivkdi cos O ) * Re A , 



(40) .> 



which is in agreement with a formula given by 
Feinberg [1946]. It is evident that the effect is very 
small when d x is greater than a few wavelengths. 
(When propagating from land toward the sea Re A 
may be replaced by (eow/O* where e =8.85X 10~ 12 
and (r g is the conductivity of the land.) 

6. Magnetic Field Near the Coastline 

In the analysis in this paper, the mutual impedance 
ratio Az m /z m is valid when the antennas at A and B 
are vertical electric dipoles. Actually, if neither A 
or B is near the boundary, the results are still appli- 
cable for other antennas which transmit or receive 
vertically polarized waves. However, if one of the 
antennas, say B is close to the boundary, the mutual 
impedance ratio will depend on whether the antenna 
is a vertical electric dipole (e.g., a whip) or a hori- 
zontal magnetic dipole (e.g., a vertical loop). For 
the general case of a loop antenna at B and a vertical 
dipole at A, it is possible to use the same method 

outlined in this paper. However, a complexity 

-> 
arises since the tangential magnetic field vector H bt 
has additional terms proportional to 1/i? 3 and these 
tend to complicate the integrations. However, if the 
direction of propagation is perpendicular to the coast- 
line, the final results may be obtained in a relatively 
straightforward manner. 

If the field incident on the boundary, at x=0, has 
the form 

E =e- ihx , 

then the resultant vertical electric field may be 
written 

E z (z)=e-**[1+F(x)], at 2=0, (41) 

where 

F(x)=^> 



6?? 



is the mutual impedance ratio given by (13) with 
C x = \ and di=x. 

The magnetic field which lias only a y component 
is obtained from the Maxwellian equation, written 
in integral form, 



H y (x)—H v (x )=ie oi I E z (x')dx' 



Jxo 



(42) 



where x is some reference distance. Now, if 
k£o>!>l, it is safe to assume that 



VoH t (x Q ) a -E.&) = -e-^o[l+FCx )]. 
Then, without difficulty one finds 



where 



G(x) 



—VoH//-- 



*[! + <?(*)], 



(43) 



(44) 



--F(x )e,- m *o-*'-{-ike mx [ x °F(x')e- ikx '<lx'. (« 
Jx 



5) 



Since F(x) is known, the function G(x) can be found 
from an integration which can be best done numer- 
ically. It can be seen that, when fe>>l, G{x)^ 
F{x) , as it must. 

It is interesting to note that, for an abrupt bound- 
ary, G(x) does not have the singular behavior char- 
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acteristic of F(x), at the boundary. While F(x) 
behaves as log x, the function G(x) behaves as 
x(\ogx)— x when |A:a?|<C!<C 1- 

7. Appendix 

For convenience in auxiliary calculations, a 
graphical plot of the functions ^(x) and g 2 (x) is 
presented in figures 7 to 10. In complex form these 
quantities are defined by 

( h ( x ) = (l/2) exp fo + ;(37r/4)]J7o (2) (!x|) for X ^0, (46) 
ft(x) = (x/2) expfe+i(W4)][Si 2) (x)-^i 2) (x)], (47) 
for x>0> an( l 

2 (x) = (x/2) exp [ix+%W*)] [Ho 2) (-x)+^i (2) (-x)L 

(48) 
for x <0. 

The mutual impedance increment Az m may then 
be calculated for any value of O by noting that 



5 =Ao[cipi(«i)+^ft(ai)]' 

where d = cos O . 



A : „ 



(4!)) 
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Figure 9. Argand plot of g 2 (x) for x>t). 
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